Spin-orbit coupling and Berry phase with ultracold atoms in 2D optical lattices 
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We show how spin-orbit coupling and Berry phase can appear in two-dimensional optical lattices 
by coupling atoms' internal degrees of freedom to radiation. The Rashba Hamiltonian, a standard 
description of spin-orbit coupling for two-dimensional electrons, is obtained for the atoms under 
certain circumstances. We discuss the possibility of observing associated phenomena, such as the 
anomalous Hall and spin Hall effects, with cold atoms in optical lattices. 
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In the last decade many experiments have studied the 
quantum motion of ultracold atoms in periodic poten- 
tials created by standing waves of light [j] . Phenomena, 
such as Bloch oscillations and Wannier-Stark ladders, im- 
possible to observe for electrons in metals due to short 
relaxation times, have been observed. Atoms in different 
internal (spin) states can interact with laser light in dif- 
ferent ways depending on its polarization. This effect has 
recently been applied to the experimental study of quan- 
tum transport of atoms in one-dimensional optical lat- 
tices in the localized regime ^ j3j , for which a wide range 
of phenomena has been theoretically proposed Q. Here 
we consider ultracold atoms dynamics in two-dimensional 
(2D) optical lattices in the itinerant regime where motion 
through the lattice is of concern. 

An intrinsically new aspect of a higher dimensional 
lattice is the possible appearance of a geometric phase 
in crystal momentum space 0, 1^ • This can arise when- 
ever time- reversal and/or spatial inversion symmetries 
are broken in the lattice, or in the presence of spin-orbit 
coupling. In solid state systems, such as ferromagnetic 
crystals for example, the geometric phase is responsible 
for the anomalous Hall effect (AHE) 7], that is the gen- 
eration of a transverse current by an electric field even 
in the absence of a magnetic field. A closely related 
phenomenon, recently proposed for spintronics applica- 
tions 8], is the spin Hall effect (SHE) [9|, that is the pro- 
duction of a transverse spin current by an electric field in 
the presence of a significant spin-orbit coupling. The ob- 
servation of the spin Hall effect with present technology 
is however still an open challenge in solid state systems. 

In this Letter, we show how spin-orbit coupling and 
geometric phase may be produced for atoms in 2D opti- 
cal lattices, and we propose experiments to explore their 
consequences. Optical lattices provide great flexibility 
with which potentials can be created and atomic quan- 
tum states prepared. By choosing the polarization of 
the beams appropriately, the internal degrees of freedom 



of the atom can be coupled to their momenta as in the 
(relativistic) spin-orbit effect for electrons in solids. In 
particular, we show how a Hamiltonian similar to the 
Rashba Hamiltonian for electrons in 2D semiconductor 
systems 10] can arise in the description of atoms propa- 
gating in an optical lattice produced by the interference 
of suitably polarized laser beams. A constant force field, 
analogous to an electric field for electrons can be pro- 
duced by accelerating the lattice Q. By studying the 
transport of the atoms in these spin-dependent lattices, 
one can observe effects similar to AHE and SHE. 

Let us consider the Hamiltonian for an atom inter- 
acting with a configuration of laser beams producing an 
electric field E(r). In general, if the detuning of the 
light frequency lo from the resonance frequency is large 
with respect to the radiative width of the excited states, 
spontaneous emission is suppressed and we can adiabat- 
ically eliminate the excited states by writing an effective 
Hamiltonian which involves only the ground state sub- 
levels [Ulli: 
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For a ground state of total angular momentum F, the 
indexes a and f3 run over the 2F -\- 1 Zeeman sublevels. 

= {e : j\d\g : f3) is the dipole matrix element be- 
tween the ground state sublevel f3 and the excited state 
sublevel 7 (of energy hio^). 

Given the analogy with electrons, in the following we 
shall focus on the case of atoms with F = 1/2. This case 
can be effectively realized using alkali atoms, cooling and 
manipulation of which have gone a long way. In particu- 
lar, ^Li atoms have a ground state hyperfine component 
of F = 1/2. Similar features are expected to occur for 
higher values of the atomic spin F as well, but we leave 
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their detailed analysis for the future. 

In the F = 1/2 case, the external potential VQ^(r) 
can be simply written in terms of a fictitious magnetic 
field B(r) coupled to the total atomic angular momentum 
operator F: 



Q/3- 
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The scalar potential V{r) is proportional to the local light 
intensity, while the vectorial field B(r) is proportional to 
the local electromagnetic spin: 



V{r) = feoE*(r).E(r) 
B(r) = -i6iE*(r) x E(r). 



(3) 
(4) 



The proportionality coefficients &o,i depend on the de- 
tails of the atomic structure as well as on the light fre- 
quency. As discussed in |Q| an effective coupling to the 
fictitious magnetic field requires the detuning from the 
excited state to be smaller than the fine structure of the 
excited state. Some consequences of the fictitious mag- 
netic field have already been investigated from many dif- 
ferent points of view, for instance Sysiphus cooling in res- 
onant optical lattices , or NMR type experiments ■ 
The fictitious magnetic field was crucial for the observa- 
tion of mesoscopic quantum tunneling in one-dimensional 
spin-dependent optical lattices 0. 

We consider a 2D optical lattice created with three 
laser beams propagating in a plane with equal angles 
between them (see Fig. [T]), a geometry already studied 
in the dissipative regime [l5| . As done in 0|, we how- 
ever restrict our attention on the dissipationless case. We 
choose the polarization of the electric fields in such a way 
that the components respectively parallel and perpendic- 
ular to the plane are the same for all three beams 
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the Qi are the wave vectors of the light beams, coefficients 
a and f3 are complex and z is the unit vector perpen- 
dicular to the plane. For this configuration, the scalar 
potential © and the fictitious magnetic field are 

V{r) = ^Focosk,r, (6) 

i 

B(r) ~ (Viz sinker + V2k^ sinker + V^z x k^ cosk^r), 

i 

where ki = q2 — qs, e^c, and the amphtudes Vi are 



K)-6o(|a|' + |/3|'cos(2V3)), 

Vi = 2bi |/3|^sin(27r/3), 

V2 = -4fei sin(7r/3)Re {a* (3) , 

V3 =26iIm(a*/3). 



(7) 



We will use M, (qV^) ^ and 6Er as the basic units of 
mass, length and energy, where M is the mass of the 
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FIG. 1: Spatial profile of the different terms of the poten- 
tial JSJ: (a) Vo term, (b) Vi term, (c) V2 term, (d) V3 term. 
The amplitude of a vector component is proportional to the 
length of the arrow. The inset in panel (a) shows the con- 
figuration of the laser wave vectors (red arrows) and the 
lattice vectors ki (blue arrows). Red dots indicate wells A, 
green triangles indicate wells B. Dashed lines show a unit cell. 

atom, qi is the light wave vector, and = fi^q^ /2M is 
the recoil energy. Experimentally, potentials with depths 
up to 20Er are achievable in the dissipationless regime. 

In what follows we are interested in the dynamics of 
an atomic distribution prepared in the lowest band of 
the potential. The different terms in the Hamiltonian 
proportional to the Vi are shown in Fig. ^ We consider 
the case in which |Vb| ■> \Vi\ , \V2\ , l^al (|a| > To- 
gether with a direct numerical calculation of the band 
structure and Berry curvature (see Fig. \^ , we consider a 
tight-binding model that allows us to obtain an effective 
Hamiltonian close to the F point. The scalar part of the 
potential creates a honeycomb lattice which is perturbed 
by the vector part. Introduction of the term with Vi 
makes the lattice asymmetric for each spin component, 
similar to the one recently studied by us ^2 and 

V3 lead to hopping processes with simultaneous spin flip- 
ping. The tight-binding approach results in the effective 
Hamiltonian 
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Since lattices for different spins are mirror images of each 
other, the on-site energies for different spins are such that 
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FIG. 2: (a) Dispersion of the lowest four bands in the po- 
tential 0. The green dashed line corresponds to the case 
in which only a scalar potential is present, Vo = 1- The blue 
solid line is for the case Vo = 1, Vi = V3 = 0.1 and V2 = [T9l| . 
(b) Berry curvature for the lowest band in the presence of an 
external field ho — 0.005. The band structure in the presence 
of ho is essentially the same as shown in (a), except for a 
splitting by « 0.01 at the degeneracy points F, K and K'. 



-Cj^ — e]j. OfF-diagonal matrix elements are 



-t (1 + e 
f(e-'^/3 + e 



i7r/3g — ikR,i 



(9) 



-ikR2 



where t is proportional to the overlap of the wave func- 
tions in different sites with the same spin, t is propor- 
tional to the overlap of wave functions at different sites 
with different spin, and are the lattice primitive vec- 
tors (see Fig. The coefficient t is determined by the 
scalar potential, while t is governed mostly by V3 terms 
(to the lowest order, terms with V2 do not affect the dy- 
namics, since the corresponding B field vanishes at the 
place where the overlap between on-site wave functions 
is the largest, see Fig.^. The parameter ho describes an 
additional external field. 

For vanishing fictitious magnetic field B and external 
field ho, the bands have a two- fold spin degeneracy at all 
points of the Brillouin zone; at points K and K', where 
bands cross, the degeneracy is four-fold. As the Vi term 
is added, the bands keep the two-fold spin degeneracy, 
but the four-fold one at K and K' is lifted and a gap is 
correspondingly opened. Spin degeneracy is then lifted 
by the inclusion of the V3 term. The band degeneracies 
at the high symmetry points can be understood from the 
tight-binding model simply by looking at the properties 
of the off-diagonal matrix elements in ©. For instance, 
at the r point one has vi^2 = and only v ^ 0, there- 
fore the upper and lower bands are two-fold degenerate 
there. An effective Hamiltonian close to the F point can 
therefore be obtained simply by treating vi and V2 as a 
perturbation. 

The spectrum of the unperturbed Hamiltonian in the 
absence of the external field consists of two doubly de- 
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For small k around the F point, the dispersion is 
quadratic in k. The corresponding two lowest eigenvec- 
tors are of the form: 



\,-) =a\AJ)+b\B,^), 
I,-) ^b\A,i)+a\B,i). 



(11) 



In the subspace spanned by these two eigenvectors the 
effective Hamiltonian once vi and V2 are taken into ac- 
count has the form: 



iJr(k) = e_/ + hoff;: - j^kya^ + k^ay 



(12) 



where cr^ are Pauli matrices and / is the unit matrix. This 
reduces to the standard form of the Rashba Hamiltonian 
if one makes a global spin-rotation about the ay axis to 
flip the signs of ax and a^ . For the model here discussed, 
the value of the spin-orbit coupling parameter 7 ~ 0.08 
can be extracted from the relative slope of the two low- 
est eigenstates. A Hamiltonian of such a form has been 
recently predicted to give both AHE "T] and SHE ^ in 
solid state systems. 

The simple analytic form of the Hamiltonian in H12() 
gives a simple expression for the Berry curvature which 
is responsible for AHE 



{hl + j^k^) 
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(13) 



The results of the Berry curvature calculation for the 
continuous potential in Eq. ((HJ are shown in Fig. |2] For 
ho — the Berry curvature is zero everywhere and sin- 
gular at the points where the bands touch. Any finite 
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FIG. 3: (a) AHE. Trajectory of the wave packet center for 
a force f in the positive y direction. Red dots are period- 
averaged positions. The initial position is indicated by the 
blue dot. (b) SHE. Distribution of the spin vertical compo- 
nent as a function of momentum when finite force f = fxX is 
applied. Red is positive, blue is negative. The blue arrows 
show effective momentum dependent field Ak close to the F 
point. 
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value of the external field ho completely removes the de- 
generacy of the bands and makes the Berry curvature to 
be peaked around the F and K' points with a finite max- 
imum and width. An effect similar to AHE can therefore 
be observed with cold atoms: a wave packet initially pre- 
pared in the lowest band with a small quasimomentum 
spread Afc <C 1 0| around the F point and accelerated 
in the FK direction performs Bloch oscillations in the di- 
rection of the drive and at the same time it drifts along 
the perpendicular direction because of the geometrical 
phase accumulated In Fig.|21(a), we show the trajec- 
tory of the wave packet center calculated with semiclassi- 
cal equations including the Berry curvature effects 0, . 
The external force has been taken as / = 0.001. In order 
to maximise the drift in a given time, the largest force 
that preserves adiabatic evolution in the first band has 
to be chosen in an actual experiment. The magnitude 
of the acceleration creating the force f is ^ 500 m/s^, a 
value already demonstrated in the framework of optical 
lattices 0. The chosen value of external field ho corre- 
sponds to ^ 1 mG. To observe a wave packet drift of 
~ lOA, the acceleration needs to be applied for ~ 50 ms. 

To observe an effect similar to SHE with cold atoms, a 
spin-sensitive measurement of the wave packet momen- 
tum distribution after it was exposed to an external force 
has to be performed. The quasimomentum spread has 
to be small enough, and the force applied for a short 
enough time, so that the wave packet stays within the 
region where (|12|l is valid. This region has a radius of 
the order of 0.2 around the F point. As discussed in 9] 
off-diagonal terms in Eq. H12|l may be thought of as a mo- 
mentum dependent magnetic field Ak = —j{kyX + k^y). 
When ho = and an external force is applied, the spins, 
initially parallel to the plane, acquire a non- vanishing 
z-component, i.e. "tilt vertically" . For an atom with 
quasimomentum {kx,ky) and a force along the x direc- 
tion f = fxX, the z-component of the spin is: 



fxky 

jk^ 



(14) 



The atoms with a positive ky thus acquire a positive z- 
component of the spin, while the atoms with a negative 
ky acquire a negative one. For the momentum spread, 
we have taken Aky — 0.1, which corresponds to the 
minimum of the first band when ho = 0, and for the 
force / = 10~^ 50m/s^). For these values, we predict 
n-z,k = 0.125. To measure the spin tilting in the vertical 
direction one should measure the momentum distribution 
of the different vertical spin components after the force 
has been applied. Spin dynamics can be frozen by sud- 
denly switching on an external ho S> ^A.k. The lattice 
potential is then adiabatically removed: this transforms 
the quasimomentum distribution into a true momentum 
one so that SHE can be observed as the motion of the 
different spin z-components in opposite directions along 
the y axis. 



In conclusion, we have considered continuing the par- 
allels between the dynamics of electrons in solids and 
cold atoms in optical potentials. We have shown that 
the coupling of the atomic internal degrees of freedom to 
the translational ones can be described in terms of the 
same Rashba Hamiltonian as used to model spin-orbit 
coupling in solids, and we have predicted the possibility 
of observing AHE and SHE with cold atoms. With re- 
spect to solid state systems, atomic ones are expected to 
provide a much cleaner environment and a easier tunabil- 
ity of the system parameters. This will allow for a deeper 
understanding of the basic physical mechanisms that may 
be important for future spintronics applications. 
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ones trapped in the lowest band. no transverse drift coming from the band asymmetry can 
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